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Abstract 

We report experiments on the rapid contact hne motion present in the early stages of capihary 

driven spreading of drops on dry solid substrates. The spreading data fails to follow a conventional 

viscous or inertial scaling. By integrating experiments and simulations, we quantify a contact line 

friction (/Uj), which is seen to limit the speed of the rapid dynamic wetting. A scaling based on 
c/i ■ 

this contact line friction is shown to yield a universal curve for the evolution of the contact line 
^ radius as a function of time, for a range of fluid viscosities, drop sizes and surface wettabilities. 
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The interest in moving contact lines is increasing due to the need for design of fluid 
applications at small scales, since these often rely on manipulation or control of two-phase 
flow. Examples of such are microfluidic systems, sintering, printing, coating and immersion 
lithography techniques to name a few. Spontaneous spreading occurs in the deposition 
and formation of micron sized drops in biomedical applications and when rewetting the 
lubricating film covering the eye. 

A generic example of dynamic wetting is the spreading of a spherical liquid drop as it 
comes in contact with a dry solid surface. Its spreading after contact is dominated by 
different physical mechanisms at various stages in the temporal evolution. If the drop radius 
is less than its capillary length, the flow is mainly driven by the interfacial energy of the drop 
and the substrate surface energy. The contact line is formed at the intersection of the drop 
liquid-air interface and the solid substrate, where the dynamic contact angle is defined as 
the angle between the liquid-air interface and the substrate. For a moving contact line, the 
interface is typically distorted near the solid surface, giving rise to a free surface capillary 
force, which may pull the contact line forward. These forces are balanced by different rate- 

n n 

limiting processes, such as viscous dissipation |1[ and inertia |2[, which all act to reduce the 
contact line speed. 

It is well known that the classical hydrodynamic theory predicts a divergence of viscous 
stress at the contact line. Therefore it might be expected that the spreading is dominated 
by the viscous dissipation in the bulk. By regularizing the viscous dissipation, a model for 
the spreading in viscously dominated wetting is established jsj. This is often referred to as 
Tanner's law where the spreading radius (r) evolves as r ~ i?(^)To, where a is the surface 
tension coefficient, R the initial drop radius and /i the viscosity. This model, which holds 
promise if the drop evolves slowly and has a shape similar to a spherical cap, has explained 
many experiments. However there are many wetting phenomena that it does not describe, 
illustrating that there are other mechanisms influencing or dominating the spreading. 

One example is the spontaneous spreading of a water drop as it comes in contact with 
a low energy substrate. Experiments indicate here that the acceleration of liquid in the 
bulk of the drop is resisting contact line motion. An inertial spreading is found to follow 
r ~ i? * ■ \a\ (p is the density), but by makirig the substrate more hydrophobic a 

different exponent for the spreading radius was found [2] . The hydrodynamic model cannot 
fully capture wetting at high capillary numbers (given by the ratio of the viscous and surface 
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tension force) 5j, and dynamic wetting experiments of viscous (IPa s) drops 6|. In the latter 
case the spreading radius was observed to increase as the square root of time (r ~ t^/^). 

Q 

De Gennes [7| postulated that there might be another non-hydro dynamic dissipative 
contribution arising from the contact line itself. This macroscopic dissipation was defined 
by a friction factor local at the contact line, which has the same units as viscosity. Others 



12| have also discussed the importance of local non-hydrodynamic effects at the contact line, 



with different interpretations of its microscopic origin. Recently [13| a friction factor was 
estimated from the molecular kinetic theory by fitting the experimental spreading radius 
for drops with different viscosity. These experimental observations are collected at much 
later time scales than presented here, and the value for this friction factor is an order of 
magnitude larger than our numerical measurements. 

By integrating experiments and axi-symmetric simulations based on the Cahn-Hilliard 
Navier Stokes equations we estimate values for the friction factor (fif) that appears 

in the free energy formulation. Theoretically, the friction factor generates a local dissipation 
at the contact line through its boundary condition. Here, particular attention is devoted to 
the very first stage of a spontaneous spreading process that is far from equilibrium. The 
experimental data cannot be rationalized as viscous or inertial effects. The data set collapses 
for a scaling law based on the numerically measured contact line friction parameter /ij, even 
for a wide range of viscosities (l-85mPa s), different drop sizes and surface energies. These 
results indicate that local dissipation at the contact line, interpreted as a contact line friction, 
is limiting spreading. 

Both experiments and numerical axi-symmetric simulations of drop spreading have been 
performed. The simulations are based on the Cahn-Hilliard Navier-Stokes equations fl^ . 
In terms of phenomenological thermodynamics one can postulate the free energy [F) for a 
binary fluid F = J (f ^(C) + f jVCp) dn + J ((a,^ - asg)g{C) + a^g) dT. The volumetric 
(fl) free energy consists of two terms representing the bulk (^\I'(C)) and interfacial energy 
(yI^C'P), respectively. \E' = — 1)^ is a double- well function with two minima, giving 
the equilibrium values of the order parameter C, as C = —1 for gas and C = 1 liquid. 
The diffuse interface width (e) is chosen to be the same as the spatial resolution in the 
experiments e = 7.5/xm. Important to note, however, is that in [8;] e has been varied one 
order of magnitude, without any noticeable change in the results or any increase in viscous 
dissipation. 
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The surface energy of the wet substrate is asi, and the dry ((Xsg). g{C) = |(2 + 3C — C^) 
is chosen to give g{l) = 1 and g{—l) = 0, thus producing the corresponding wet or dry 
surface energy of the substrate. 

By making a variation in F with respect to the concentration, one obtains an expres- 
sion for the chemical potential {6F/6C). If accounting for the effects of convection of the 
concentration, that would equal the flux due to gradients of the chemical potential, the 
Cahn-Hilliard equation is recovered, which along with the Navier Stokes equations forms a 
theoretical basis for modeling of wetting [l^ with a no-slip on the wall. 

By retaining any perturbation in the concentration at the wall, a general wetting bound- 
ary condition for the concentration at the solid surface appears [la], 

dC 

e/"/^ = -e^VC ■ n + acos{e,)g'iC). (1) 

We interpret here /xj as a friction factor at the contact line. 6^ is the equilibrium contact 
angle. 

Experiments of spontaneously spreading drops have been carried out through high-speed 
imaging (ISOkfps) for different viscosities and coatings (oxide, silane, teflon) on Si- wafers. 
The viscosity was changed by using different glycerin-water mixtures, of glycerin mass- 
fraction 0%, 50%, 65%, 72.5%, 82.5% , corresponding to viscosities [1, 6.6, 14, 31, 85] mPa 
s, respectively. The different viscosities do not give any significant change in equilibrium 
contact angle (±2°), which were measured as 9e = [20°, 60°, 109°] for oxide, silane and teflon 
coatings. 

The axi-symmetric Cahn-Hilliard Navier Stokes simulations mimic the experiments using 
the same material properties (density, viscosity, surface tension and equilibrium angle) as 
measured from experiments. To obtain the experimentally observed spreading behavior, an 
additional dissipation at the contact line was necessary through a non-zero Hf [l^. /i/ was 
determined by obtaining a direct agreement between simulations and experiments, enabling 
a direct measurement of fif even in the presence of other contributions such as viscosity and 



inertia 



15| . The values for fif are reported in table I for all the surfaces and viscosities. 



A non-monotonicity in fif is observed for pure water for the Si02 and silane coating, the 
same dependency was reported in [14] when comparing with similar experiments [3]. We 
can at the present time not explain this non-monotonicity for pure water. Fig. 1 shows 
the excellent agreement between simulations and experiments for a water and a glycerin- 
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a) b) c) 



FIG. 1: Panel a) Illustrates the initial condition for the experiments and the numerical 
simulations, where a drop held at the tip of a needle is brought into contact with a dry 
solid substrate. The dashed box in the figure shows the field-of-view in the experiments. In 
panels b) and c) are shown the drop shape near the substrate, at times t = 0.15ms, 
t = 0.60ms and t = 0.90ms, after initial contact. Each panel shows a composite of 
experiment (left) and simulation (right). The black solid line in the right half that is 
plotted on top of the simulation result, illustrates experimental interface shape, b) A water 

drop spreading on an oxidized Si-wafer {9e = 20°, viscosity fiH20 = ImPa s). c) Glycerin 
82.5% drop spreading on an oxidized Si-wafer {9^ ~ 20°, viscosity iJigiycering,2^% = 85mPa s). 



Mass fraction glycerin 


0% 


50% 


65% 


72.5% 


82.5% 


Si02 [Pa s] 


0.15 


0.33 


0.51 


0.66 


1.02 


Silane [Pa s] 


0.17 


0.26 


0.33 


0.41 


0.80 


Teflon [Pa s] 


0.07 


0.06 


0.09 


0.10 


0.19 



TABLE I: Values for the contact line friction parameter [Pa s] for different viscosities 
and substrates {Si02, Silane, Teflon) measured from the numerics. 



water 82.5% drop with an initial radius R ~ 0.5mm. See also Supplemental Material at 



20l |. Fig. la shows the initial condition in the experiments and simulations, and the field of 



view in the experiments (dashed box). The same window was extracted from the numerics, 
however the whole drop was simulated. 

Fig. 2a shows how the radial position of the contact line evolves in time for drops with 
different initial radii and for different viscosities on the oxidized Si-wafer. The markers 
represent the mean value after several realizations of the experiments (minimum of four) 
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FIG. 2: The spreading radius in time on an oxidized Si- wafer for two drop sizes 
R ~ (0.3 ± 0.02)mm (hollow markers) and R ~ (0.5 ± 0.02)mm (filled markers) for 
different mass fractions glycerin as indicated in the legend, (a) Dimensional units, (b) 
Viscous scahng. (c) Inertial scahng. (d) Contact hue friction scaling. 



and the data set has been reduced for clarity. One observation to be made in fig. 2a is that 
the viscosity as well as the drop size influences the spreading. 

Fig 2b shows the same data, with the contact hne radius scaled with initial drop radius R 
and the time scale with a viscous capillary speed a/fi. The capillary speed a/ fi is 73m/s for 
water and 0.75 m/s for 85% glycerin-water. However, as is evident from fig 2b, this scaling 
fails to collapse the data, so the viscous contribution does not seem to be the limiting factor 
in this situation. An alternative would be an inertial scaling of time based on an inertial 
capillary velocity scale ■\/a/{pR), as shown in flg 2c. As is evident here, this scahng does not 
capture the essential dynamics either, and we conclude that neither inertia or bulk viscosity 
is the limiting factor for spreading in our experiments. 

The remaining possibility is a capillary velocity based on the contact line friction discussed 
above and quantified in table 1. A representative velocity in this case can be found either 
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from equation [2] or from dimensional analysis to be u* = cr/fif- Introducing the values for 
cr and fif from table 1 gives a speed of u^y^ ~ 4.8m/s for water and Wg^^^ ~ 0.6m/ s for 85% 
glycerin. By scaling time with R/u*, we do obtain a collapse of data, for the entire range of 
viscosities and drop sizes, see fig. 2d. The scattered dimensional plot represented in fig. 2a is 
reduced to nearly a single spreading curve. Fig. 2 shows only results for the Si02 surface, 
but similar results are also obtained for the other solid surface coatings. 

/i/ is determined by adjusting it in simulations so that the mean spreading radius agrees 
with that of several experiments performed using the drop radius 0.5mm. It should be 
noted that the adjustment of this single parameter achieves excellent agreement for the 
entire drop shape, over the whole spreading event. We have also varied the drop size in 
additional experiments, which has a significant influence on the spreading radius (see fig. 2a). 
As shown in fig 2d, the data for both drop sizes collapse excellently when using a scaling 
of time according to at/{R * jif). The value of /i/ is thus independent of drop size, and 
this indicates that it is an intrinsic material property of the surface in combination with the 
wetting liquid. 

Fig. 3a shows the non-dimensional collapse of data for the three surface coatings for dif- 
ferent drop sizes and viscosities. By representing the dimensionless curves in fig. 3a in loga- 
rithmic axis, we observe that the radii follow the same slope independent of the solid surface 
at the early stage of the partial wetting, see fig. 3b. This is indicating that the governing 
physical mechanism is indeed the same for the different solid surfaces. From fig. 3b it is 
clear that the spreading radius evolves as ~ ("^)^- ^ similar relationship is expected 
in a diffusion process, where in this context aR/fif would represent a diffusion coefficient. 
This could indicate that a diffusive process is taking place at the contact line, which was 



suggested by 



10| from rapid wetting simulations using molecular dynamics. However, a 



detailed study at the nanoscale would be needed to verify this. In the first stage of the 
spreading, for non-dimensional time < 1, the experiments cannot be fully captured by the 
hydrodynamic theory through Tanner's law r = -R(^)^ or by the molecular kinetic theory 



17, 



18| . We have for clarity inserted the slope predicted from Tanner's 



that predicts r ^ t'l 
law in fig. 3b. 

In fig. 3b a distinct transition between the 1/2 slope and a much more gradual slope (~ 
1/10) is observed around non-dimensional time 1. This might be an indication of the transi- 
tion between contact line friction dominated spreading and another slower spreading regime. 



7 



We assume here that the second regime is viscously dominated spreading given by Tanner's 



law and makes this equal to the contact line friction dominated spreading r = R^J a 
distinct transition time {tt) between the two regimes is obtained. In dimensional scales this 
becomes tt = ~~(^)^ or in non-dimensional time (r), r = = (^)^. Introducing the 
material properties in the expression for r we notice that a physically reasonable transition 
time is obtained and in very good agreement with the experimental results presented in fig. 3. 
For example, the dimensionless transition time for water and 85% glycerin is on the oxide 
surface found to be tq% = 3.5 and T8^% = 1.05, respectively. 

An analytical function can be derived for the contact line velocity {ud), based on the 
boundary condition given in eqU] if the equilibrium profile for the concentration across the 
interface is introduced and some algebra is performed jl9|. 

a cos{6e) - cos{6) 

Ucl = TT^T (2) 

where 6 is the dynamic or apparent contact angle. Eq. [2] is different from other expressions 
for the contact line velocity previously reported in the literature 7|] by that it is divided by 
sin(^) which makes the expression diverge at angles 0° and 180°. This function is assumed 
to only be valid when the local dissipation at the contact line dominates. At these extrema, 
other mechanisms such as inertia or bulk viscous friction are expected to regularize the solu- 
tion, sin (6') gives a non- negligible contribution to the function and introduces an additional 
non-linearity. 

In fig. 3b it is clear that the spreading radius evolves as a function r ~ ^(^^)^ ■ 
taking the time derivative of this expression we find that the contact line speed should be 
proportional to ~ t^^. To evaluate the analytical expression for the contact line velocity 
given in eq.(2), we use the experimental data for the dynamic contact angle for the data 
presented in fig. 2a for the different viscosities and drop sizes, as they evolve on the oxidized 
wafer. We define the dynamic contact angle between the tangent along the contoured in- 
terface (interpolated at a fixed height of 7 pixels from the wall) and the solid substrate, on 
the liquid side. The dynamic contact angle measurements are found to be fairly insensitive 
to changes in interpolation height, as long as this height is chosen to be less than the local 
radius of curvature at the contact line. 

Fig. 4 shows that the expression given in eq.2 indeed gives a slope for the contact 



line speed as Ud ~ ^- This indicates that ^"'^^t.'i^e"'^''^'' ~ R/''^, which from 
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eq.2 recovers the experimentally observed behavior presented in fig.2d and fig. 3. The 
inset in fig. 4 shows the predicted contact line speed using the linearized function from 
molecular kinetic theory umkt = {(^/l^f) " (cos(6'e) — cos(6')) Since we are interested 

in the slope for the contact line speed in time, we assume fif to be the same in umkt 
as reported in table I. One clear observation to make from the inset in fig. 4 is that at 
non-dimensional time < 2.4 the slope for the contact line speed predicted from molecular ki- 
netic theory umkt = {l^f / o')umkt does not agree with the experimental observation in fig.3. 

In summary we have shown that spreading experiments and simulations for a wide 
range of viscosities, on substrates with very different wetting properties, all exhibit a 
universal spreading behavior if contact line friction dominates the spreading. An expression 
for the contact line radius is proposed for this spreading regime as r ~ R{-^)^ . The 



analytical contact line velocity from phase field theory, where the dynamic contact angle is 
primary input, predicts the same slope for the spreading as found directly in experiments. 
A criterion is proposed to determine the dominance of contact line friction or viscosity in 
spreading, which is found in good agreement with the experiments. We hope that these 
results can help rationalize spreading phenomena that falls beyond classical hydrodynamic 
theory, and gives a phenomenological explanation for such physics. 




9 




log(at/(R*^)5 



FIG. 3: Non-dimensional spreading radius based on a contact line friction scaling on the 
different substrates (oxide (back), silane (red), teflon (blue)). Hollow markers denote 
[R = 0.3mm) and filled markers {R = 0.5mm). (a) Linear axis, (b) Logarithmic axis. 
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FIG. 4: The main figure shows the dimensionless contact hne velocity function from phase 
field theory Ud = ^Uci = (cos(^^e) — cos{9))/ sin(6'). Inset shows the velocity predicted from 
the linearized molecular kinetic theory umkt = {fJ'f /o')umkt = cos{9e) — cos{9). The input 

in these two functions are the experimentally measured dynamic contact angle 6 for two 
different drop sizes on the oxidized Si-wafer. The mass fraction glycerin is indicated in the 
legend. Hollow markers denote small {R = 0.3mm) and filled markers large drops 

{R = 0.5mm). 
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